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Abstract 

We report on Suzaku and Chandra observations of the young supernova remnant 
CTB37B, from which TeV 7-rays were detected by the H.E.S.S. Cherenkov telescope. 
The 80 ks Suzaku observation provided us with a clear image of diffuse emission and high- 
quality spectra. The spectra revealed that the diffuse emission is comprised of thermal and 
non-thermal components. The thermal component can be represented by an NEI model 
with a temperature, a pre-shock electron density and an age of 0.9±0.2 keV, 0.4±0.1 cm -3 
and 650^300° yr, respectively. This suggests that the explosion of CTB37B occurred in 
a low-density space. A non-thermal power-law component was found from the southern 
region of CTB37B. Its photon index of ~1.5 and a high roll-off energy (>15 keV) indi- 
cate efficient cosmic-ray acceleration. A comparison of this X-ray spectrum with the TeV 
7-ray spectrum leads us to conclude that the TeV 7-ray emission seems to be powered 
by either multi-zone Inverse Compton scattering or the decay of neutral pions. The point 
source resolved by Chandra near the shell is probably associated with CTB37B, because 
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of the common hydrogen column density with the diffuse thermal emission. Spectral and 
temporal characteristics suggest that this source is a new anomalous X-ray pulsar. 

Key words: acceleration of particles — ISM: individual (CTB37B) — ISM: supernova 
remnants — X-rays: ISM 

1. Introduction 

Supernova Remnants (SNRs) are one of the most promising acceleration sites of cosmic rays. 
In fact, ASCA detected synchrotron X-ray emission from the shell of SN 1006, which unambiguously 
indicates the acceleration of electrons up to ~100 TeV (Koyama et al. 1995). Following this discov- 
ery, the synchrotron X-ray emission has been discovered from a shell of a few more SNRs, such as 
RX J1713.7-3946 (Koyama et al. 1997) and RCW 86 (Bamba et al. 2000). On the other hand, TeV 
7-rays have also been detected from some non-thermal shell-type SNRs. The radiation of TeV 7-ray 
is explained by either (1) Inverse-Compton scattering (IC) of cosmic microwave background photons 
by the same high energy electron giving rise to the X-ray synchrotron emission or (2) the decay of 
neutral pions that are generated by collisions between high energy protons and dense interstellar mat- 
ter. The ratio of fluxes between the TeV 7-ray and the X-ray provides the magnetic field intensity 
as long as one assumes that the TeV 7-ray is produced through the IC mechanism. Utilizing this 
characteristic, Matsumoto et al. (2007) found that the TeV 7-ray from HESS J1616 — 508 is likely the 
result of the proton acceleration, because the non-detection of X-ray using the Suzaku XIS provides 
much weaker magnetic field than the interstellar average. 

Although the evidence of particle acceleration has accumulated rapidly, our knowledge is still 
limited on what sort of conditions are necessary for SNRs to accelerate particles. A breakthrough 
may be brought about by searching SNRs from which the TeV 7-ray emission is already detected for 
thermal emission systematically, since the thermal emission provides us with a lot of information on 
the environment such as temperature, density, and age of the plasma. 

CTB37B locates at (I, b) = (348.°7, +0.°3) with a distance of 10.2 ± 3.5 kpc (Caswell et al. 
1975). This region is one of the most active regions in our galaxy where star burst activities, a number 
of shell structures probably associated with recent SNRs (Kassim et al. 1991), and OH maser sources 
(Frail et al. 1996) are detected in radio band. TeV 7-ray emission is also detected by the H.E.S.S. 
observation (Aharonian et al. 2007). In spite of the evidence of the high activities in other wave 
bands, X-ray observations have been relatively poor. Only ASCA (Tanaka et al. 1994) has detected 
a part of CTB37B at the edge of the field of view of the Gas Imaging Spectrometer (GIS) (Ohashi 
et al. 1996; Makishima et al. 1996) in the course of the galactic plane survey (Sugizaki et al. 2001; 
Yamauchi et al. 2008). Although the statistics are limited and the response of the GIS is not qualified 
at the pointing position of CTB37B, Yamauchi et al. (2008) represents that the fit of a power law to 
the GIS spectrum results in a steep photon index of ~4.1, whereas fit of an optically thin thermal 
plasma model requires a high temperature of ~ 1.6 keV. These results strongly suggest that the X-ray 
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spectrum is a mixture of a non-thermal power law and an optically thin thermal plasma emission. In 
addition, Aharonian et al. (2008) resolved a bright point source located near the shell of CTB37B 
from the diffuse emission by Chandra, although its spectral parameters are not constrained very well 
because of short exposure time. 

In order to take an image and high quality spectra of CTB37B, we have carried out an ob- 
servation of CTB37B with Suzaku. We also refer to the Chandra archival data to include the spatial 
structure and to compare them to our Suzaku data. In § 2, we present the observation log and data re- 
duction method. Image analysis is presented in § 3. Spectrum analysis and timing analysis are shown 
in § 4 and § 5, respectively. We have really detected both the thermal and non-thermal power-law 
components from CTB37B. Discussions are made in § 6 on the nature of the thermal and non-thermal 
component as well as the point source. Finally we summarize our results in § 7. 

2. Observation and Data Reduction 

2.1. Suzaku Observation 

CTB37B was observed with Suzaku (Mitsuda et al. 2007) during 2006 August 27-29. The 
nominal pointing position was (RA, Dec) = (17 h 13 m 57 s , -38°12'15", J2000). Suzaku is equipped 
with two kinds of X-ray detectors; one is the Hard X-ray Detector (HXD; Takahashi et al. 2007; 
Kokubun et al. 2007), which is a non-imaging type detector and is sensitive in the 10-600 keV band. 
The other is the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007), which is an X-ray CCD 
camera mounted on the focal plane of the X-Ray Telescope (XRT; Serlemitsos et al. 2007). In to- 
tal, there are four modules of the XIS, three of which are Front-Illuminated (FI) CCDs, which are 
hereafter referred to as XIS-0, 2, and 3, and the other one is a Back-Illuminated (BI) CCD, which is 
referred to as XIS-1 . The XRT has a point-spread function (PSF) of a Maltese-cross shape with a core 
radius of ~15" accompanied by an outskirt extending a few arcmin. The half-power diameter (HPD) 
of each telescope is ~2 arcmin. We concentrate on the XIS data in this paper because the HXD has 
no imaging capability and hence there remains a large systematic error in estimating the flux from 
CTB37B. 

The XIS was operated in the normal full-frame clocking mode with neither burst nor window 
options and SCI-off. The editing mode was 3 x 3 in low and medium data rates and 5 x 5 in 
high and super-high data rates. In analysis, we employed the data processed with the revision 1.2 
pipeline software, and used the HEADAS software (version 6.2) and XSPEC (version 11.3.2) for 
the data reduction and spectral analysis, respectively. We applied the charge-transfer inefficiency 
(CTI) correction by ourselves with the xispi software and CTI parameters of 2006-08-23. After the 
screening of the data, the effective exposure time of 80 ksec in total. The response matrix files (RMF) 
and ancillary response files (ARF) were made using xisrmfgen and xissimarfgen (Ishisaki et al. 2007) 
version 2007-09-22 under the assumption that the emissions are from point source. 
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2.2. Chandra Observation 

Chandra observation was performed on the 2 th February 2007 with the Advanced CCD 
Imaging Spectrometer (ACIS). Chips 10, II, 12, 13 S2 and S3 were used. The angular resolution 
is ~0.5 arcsec which correspond to the CCD pixel size. The data reduction and analysis were made 
using the Chandra Interactive Analysis of Observations (CIAO version 3.4, CALDB version 3.3.0). 
The total exposure time is 26 ksec after screening the data. 

3. Image Analysis 

3.1. Suzaku Images 

Fig. 1 shows Suzaku XIS images in 0.3-3.0 keV and 3.0-10.0 keV. They are created by com- 
bining those from all the four XIS modules and smoothed with a Gaussian with a = 12 arcsec, which 
is close to the XRT core size and effective in highlighting the diffuse emission. The source that locates 
at (/, b) ~ (348.°68, 0.°37) appears as the brightest source both in the soft and hard bands. Another 
source extending to the south of the brightest source, at (I, b) ~ (348.°63, 0.°32) seems to be a diffuse 
source and manifests itself only in the band above 3 keV. In addition to these sources, a point source 
is detected at (I, b) ~ (348.°56, 0°33) in the band below 3 keV. The sky position is consistent with that 
of the point source 1RXS J171354.4-381740 listed in the ROSAT Bright Star Catalogue (Voges et al. 
1999). In order to investigate these sources separately, we defined the following photon-integration 
regions (see Fig. 1) for the spectral analysis. Region 1 is the green circle with a radius of 2.'6 centered 
at the intensity peak of the brightest source. Region 2 is the blue ellipse with a major and minor axis 
of 2.'5 and lfl, respectively, which is centered at the second diffuse source. Region 3 is the circle 
colored in magenta with a radius of 1.'3. The other three regions with the same colors but with dashed 
lines define those collecting the background events. We set these background regions by taking into 
account the telescope vignetting. 

3.2. Chandra Images 

Fig. 1(c) shows the Chandra image in the 0.3-10.0 keV band corrected for the telescope vi- 
gnetting and smoothed with a Gaussian with a = 6 arcsec to see the point sources clearly. In total, 
18 sources are detected above 5a confidence level (Aharonian et al. 2008). The panels (d) and (e) 
are images in 0.3-3.0 keV and 3.0-10.0 keV, respectively, after removing the point sources. They 
are smoothed with a Gaussian with a a = 40 arcsec to highlight the diffuse emission. The definition 
of the colored regions is same as those drawn in the Suzaku images. The brightest source in region 
1 is a point source, which is defined as source A, locating at (I, b) = (348°681, 0.°371). Although 
most of the flux from region 1 originates from the point source, it is revealed from Fig. 1(d) that 
region 1 is accompanied by a diffuse emission. These two components are mixed in the Suzaku im- 
ages. As expected from the Suzaku images in Fig. 1(a) and (b), region 2 is brighter than region 1 in 
the band above 3 keV. In region 3, there is the second bright source (defined as source B, which is 
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Fig. 1. Images of CTB37B in the galactic coordinates. Panels (a) and (b) are Suzaku images in 0.3-3.0 keV and 
3.0-10.0 keV, respectively, which are smoothed with a Gaussian with a = 12 arcsec. Panel (c) is the Chandra 
image in the 0.3-10.0 keV band being smoothed with a Gaussian with a = 6 arcsec. Panels (d) and (e) are the 
Chandra images in 0.3-3.0 keV and 3.0-10.0 keV, respectively. After removing point sources, we smoothed them 
with a Gaussian with a a of 40 arcsec. Solid circles in green, blue, and magenta are the integration region of source 
photons which are named as region 1 through 3 in this order. Region 1 is a circle with a radius of 2. '6, region 2 is an 
ellipse with a size of l.'l x 2 '5, and region 3 is a circle with a radius of 1 '3. The dashed regions are corresponding 
background-integration regions. 



Table 1. Count rates of sources and diffuse emission of Chandra data. 



Region Energy band Point source count rate * Diffuse emission count rate 









[10 2 counts s 1 ] 


[10 2 counts s 1 ] 


res 


;ion 1 


0.3-3.0 keV 


2.7 ±0.1 


3.3 ±0.2 






3.0-10.0 keV 


2.4 ±0.1 


0.64 ±0.26 






total 


5.0 ±0.1 


3.9 ±0.4 


res 


;ion 2 


0.3-3.0 keV 




0.43 ±0.15 






3.0-10.0 keV 




0.93 ±0.19 






total 




1.4 ±0.2 


res 


;ion 3 


0.3-3.0 keV 


2.3 ±0.1 








3.0-10.0 keV 


0.024 ±0.021 








total 


2.3 ±0.1 





Note — All errors are at Icr confidence level. 



1RXS J171354.4-381740) locating at (/, b) = (348.°561, 0.°332), and no diffuse emission is associated 
in 0.3-3.0 keV and 3.0-10.0 keV band. Table 1 summarizes the background subtracted source count 
rates separately for point and diffuse sources in region 1, region 2 and region 3. 

3.3. Correlation with Other Energy Band 

Fig. 2 shows brightness contours of radio at 1.4 GHz in blue and of TeV 7-ray with H.E.S.S. 
in green, overlaid on the gray scale image of the Suzaku in the 0.3-10.0 keV band. Source A and B 
resolved by Chandra are represented by the filled red boxes. The radio image is taken from the NRAO 




20 40 60 
(cnt) 



Fig. 2. Gray scale Suzaku X-ray image with radio and TeV 7-ray contours (blue and green, respectively) overlaid. 
The blue radio contours are in logarithmic scale from 0.016 to 1.3 Jy/beam. The green TeV contours are in linear 
scale from 30 to 50 count/0.9 arcmin 2 stepped by 4. Red boxes are source A and B, and green cross is the center of 
TeV 7-ray emission. 
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VLA Sky Survey (NVSS) database 1 (Condon et al. 1998). The X-ray emission well conforms with 
the shell in radio. Particularly, the diffuse X-ray source detected by Suzaku in region 2 is associated 
with the southern radio sub-peak. On the other hand, the peak of TeV 7-ray emission and source A 
seems to be separated from the radio shell. Note, however, that, due to limited spatial resolution of 
H.E.S.S, the apparent TeV source morphology is consistent with that of the radio shell whose radius 
is 4!5 (Aharonian et al. 2008). 

4. Spectral Analysis 

In this section, we present results of spectral analysis of the three regions described above. We 
adopt the metal composition of Anders & Grevesse (1989) as the solar abundance. Spectral fits are 
carried out with XSPEC. We always adopt an ancillary response file (ARF) for a point source, since 
the sizes of region 1 and 2 are so small that the resultant spectral parameters including the flux will 
differ only by ~1% from the case if we take into account the spatial extent. The errors quoted are 
always at the 90% confidence level. 

4.1. Region 1 

The Chandra spectrum of source A is shown in Fig. 3(a). In extracting this spectrum, we 
took a circular integration region with a radius of 3" centered at the source. We made no background 
subtraction. Since there is no apparent emission lines, we attempted to fit a power-law model un- 
dergoing photoelectric absorption ("phabs" model in XSPEC) to the spectrum. The best-fit model is 
overlaid in the upper panel of Fig. 3(a) as the histogram. The best-fit parameters are summarized in 
table 2. The photon index, the hydrogen column density, and the intrinsic flux in the 2.0-10.0 keV 
band are 3.2±g;|, 4.0 (±0.6) x 10 22 cm 2 , and 1.8 (±0.2) x 10~ 12 ergs chi^V 1 , respectively, which 
are consistent with those of Aharonian et al. (2008). 

Fig. 3(b) is the background-subtracted spectrum of Suzaku region 1. The black and red 
crosses represent the data points from the sum of the FI CCDs and those of the BI CCD, respec- 
tively. Although there is no sign of Fe Ka line in the 6-7 keV band, we have obviously detected 
Ka emission lines from He-like Mg (1.34 keV), Si (1.86 keV), and S (2.46 keV) as well as a K0 
emission line from He-like Si (2.18 keV). This means that the spectrum includes a thermal emission 
component. We thus tried to fit the Suzaku spectrum with a model composed of a power law repre- 
senting source A and a non-equilibrium collisional ionization plasma emission model ("vnei" model 
in XSPEC; Borkowski et al. 2001; Hamilton et al. 1983; Borkowski et al. 1994; Liedahl et al. 1995) 
undergoing photoelectric absorption with a common AT H . In the fitting, we set abundances of Mg, Si, 
and S free to vary. The other abundances are fixed as solar abundance. The best-fit parameters are 
summarized in Table 2, and the best-fit models are displayed in Fig. 3(b). The fact that no iron Ka 
emission line was detected can be attributed that the non-thermal component dominates the spectrum 
in the energy band above 3 keV. The reduced \ 2 of 1 .06 implies that the fit is acceptable at the 90% 

1 http://www.cv.nrao.edu/nvss/ 
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confidence level. The temperature and the ionization parameter of the "vnei" component are obtained 

to be kT = 0.89io;?? keV and n e t [cm~ 3 s] = 3.5±i 3 ! x 10 10 , respectively. On the other hand, the photon 



3.0 ± 0.2 and the intrinsic flux is 3.3±° i 



x 10 12 ergs cm 2 s x . 



index of the power-law model is V 
The photon index is consistent between Suzaku and Chandra, whereas the flux with Suzaku seems 
to be greater than that with Chandra by a factor of ~1.8. We remark that, since there is no sign of 
emission lines in the Chandra spectrum of source A, the thermal component should be extended. 



Chandra source#l 




phabs*(vnei+power) 
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Fig. 3. (a) Chandra spectrum integrated from a circular region with a radius of 3" centered on source A. The best-fit 
model represented by a power law (histogram) is overlaid, (b) Suzaku spectra of region 1 from the FI-CCDs (black) 
and the BI-CCD (red) with a model composed of a "vnei" and a power law model. The best-fit parameters are 
summarized in table 2. 



4.2. Region 2 

The Suzaku spectra of region 2 together with best-fit model and residual are shown in Fig. 4. 
The basic features such as the He-like Si Kct emission line and no apparent sign of Fe emission line are 

phabs*(vnei+pwgpwrlw+pwglwrlw) point source fix 
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Fig. 4. The XIS spectra of region 2 with a model composed of a "vnei" and two power-law components, of which 
one represents contamination from source A. The normalization of "vnei" model and all parameters of the other 
power law are set free to vary. 
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Table 2. Best-fit parameters of the region 1 spectra 



Parameters Chandra source A Suzaku 

Power Law 
Photon Index 
Intrinsic Flux a 
VNEI 

Temperature [keV] 
abundance b Mg 
Si 
S 

n e t c 

E.M. d 
7V H e 
X 2 /d.o.f 

a Flux in the 2.0-10.0 keV band in the unit of 10~ 12 ergs cm" 2 s _1 . 
b Abundance ratio relative to the solar value (Anders & Grevesse, 1989). 
c Ionization time-scale in the unit of 10 10 s cm~ 3 , where n e and t are the 
electron density and age of the plasma. 

d Emission measure E.M. = J n e njjdV ~ n 2 V in the unit of 10 58 cm -3 , 
where n e and V are the electron density and the plasma volume. The 
distance to CTB37B is assumed to be 10.2 kpc (Caswell et al. 1975) 
e Absorption column density in the unit of 10 22 cm~ 2 . 

similar to those in region 1 . In addition to the thermal and non-thermal components, we need to take 
into account possible contamination from source A which is brightest in region 1 . In fitting the region 
2 spectra, we thus first tried a "vnei + power law (1) + power law (2)" model, where the power law (1) 
accounts for the source A contamination, and the power law (2) represents a non-thermal component 
dominating the high energy band image in region 2 (Fig. 1). The vnei component, on the other 
hand, represents the contamination from region 1, and intrinsic thermal emission from region 2 if any. 
Note that the region 2 spectra are statistically poorer than those of region 1 . Hence we have fixed the 
temperature, the abundances of Mg, Si, and S, the ionization parameter of the vnei component and the 
photon index of power law (1) at the best-fit values obtained in the region 1 fit, which are summarized 
in table 2. The flux between 2.0-10.0 keV of power law (1) is fixed at 3.3 x 10~ 14 ergs cm" 2 s" 1 , which 
is ~1% of the source A flux, on the basis of the vignetting function (Serlemitsos et al. 2007). As a 
result, the free parameters are the hydrogen column density, the normalization of the vnei component, 
and all parameters of power law (2). The result of the fit is shown in Fig. 4, and the best-fit parameters 
are summarized in table 3. The E.M. of the vnei component is ~ 1 1 (4-20) % of the region 1 best-fit 
value. On the other hand, we simulated the contamination of the vnei component from region 1 using 



3.2^3 3.0 ±0.2 

1.8 ±0.2 3.3±g:| 



89+ 21 

0.40±°;£ 
1.0±0.6 

2 1+ 1 ' 6 
4.0 ±0.6 3.6t°;2 
14.6/18 176.1/166 
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smoothed Chandra image, and found the contamination is ~ 7%. The thermal component apparent in 
region 2 spectra can therefore be entirely regarded as the contamination from region 1 , and the upper 
limit of the E.M. intrinsic to region 2 is 13% of that of region 1, or 2.7xl0 57 cm" 3 . The photon 
index of power law (2) results in Y = 1.5 (±0.4) with the reduced x 2 of 0.36. It is remarkable that the 
X-ray photon index T = 1.5 is consistent with the standard radio energy index of non-thermal SNRs 
a = 0.5. 

We next replaced the power-law (2) component by an "srcut" model, which simulates a syn- 
chrotron spectrum from an exponentially cut off power-law distribution of electrons in a homogeneous 
magnetic field (Reynolds 1998; Reynolds & Keohane 1999). According to the Green's catalogue 2 , the 
radio spectral index (a) is 0.3 with a flux at 1 GHz of 26 Jy. This small index, however, is probably 
due to contamination of thermal emission. A similar situation has been reported for 30 Dor C (Smith 
& Wang 2004). We thus fixed a at 0.5, which is the typical value of the SNRs in the radio band, and 
set the flux at 1 GHz free to vary. As a result, the normalization of the srcut model is obtained to be 
1 .43 mJy at 1 GHz. This is much smaller than the radio flux 26 Jy at 1 GHz. Note, however, that this 
radio normalization is an integration of the entire radio image of CTB37B, part of which is, however, 
spilt out of regions 1 and 2. Moreover, the radio flux encompassed in region 2 is only ^2% of that in 
region 1. We therefore estimate the radio flux within region 2 to be ~1% of the total, or « 300 mJy at 

1 GHz. Even after this correction, simple extrapolation of the srcut model well fit to the X-ray spectra 
to the radio band is much smaller than the observed radio flux. We guess that the flux in the radio 
band is dominated by thermal emission. The resultant normalization (emission measure) of the vnei 
component does not change within the error. The reduced \ 2 value is nearly the same as that of the 
power-law fit. The lower limit of the roll-off energy is obtained to be 14.8 keV. We confirmed that the 
roll-off energy does not change drastically if we varied a in the range 0.3-0.7. Based on these results, 
it is possible to interpret that the spectrum of the non-thermal component extends from X-ray to radio 
with an energy index of 0.5 in the radio band. 

4.3. Region 3 

Fig. 5 shows the background-subtracted spectra of region 3. As indicated by the images in 
Fig. 1(a), X-ray flux is detected only below ~3 keV. Since the absorption is apparently weak and 
there is Fe-L hump in the 0.7-0.9 keV band, this source seems to be a foreground point source, 
probably an active star. We thus adopted a model composed of a thin thermal collisional equilibrium 
plasma emission model ("mekal" model in XSPEC; Mewe et al. 1985; Mewe et al. 1986; Liedahl 
et al. 1995; Kaastra et al. 1996) multiplied by photoelectric absorption, and fitted this model to the 
spectra in the 0.5-2.0 keV band. The result is shown in Fig. 5, and the best-fit parameters are listed 
in Table 4. Note that the fit residuals exhibit different behavior in the 0.7-1.0 keV band between the 
FI and BI CCDs. This is probably attributed to calibration uncertainty. 

2 http://www.mrao.cam.ac.uk/surveys/snrs/ 
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Table 3. Best-fit parameters of the region 2 spectrum 



Parameters 



VNEI + Powerlaw + Powerlaw VNEI + Powerlaw + srcut 



VNEI 

Temperature [keV] 
abundance a Mg 
Si 
S 

n e t h 

EM. c 
Power Law (1) 

Photon Index 

Intrinsic Flux d 
Power Law (2) 

Photon Index 

Intrinsic Flux d 
srcut 

alpha 

roll-off E [keV] 
Normalization d 
Intrinsic Flux e 

phabs N H f 

X 2 /d.o.f 



0.89 (fix) 
0.61 (fix) 
0.40 (fix) 
1.0 (fix) 
3.5 (fix) 



0.23 



+0.18 
-0.15 



3.0 (fix) 
0.033 (fix) 

1.5 ±0.4 

no+0.07 

U - /8 -0.08 



o c+0.5 
J - J -0.7 

17.5/48 



0.89 (fix) 
0.61 (fix) 
0.40 (fix) 
1.0 (fix) 
3.5 (fix) 



0.23 



+0.14 
-0.15 



3.0 (fix) 
0.033 (fix) 



0.5 (fix) 
>14.8 
1.4 
0.78 

o c+0.5 
J - J -0.7 

17.5/48 



3 Abundance ratio relative to the solar value (Anders & (irevesse, 1989). 

b Ionization time-scale in the unit of 10 10 s cm -3 , where n e t and ; are the electron density and age of the plasma. 
c Emission measure E.M. = J n e nndV ~ n^V in the unit of 10 58 cm -3 , where n e and V are the electron 
density and the plasma volume. The distance to CTB37B is assumed to be 10.2 kpc (Caswell et al. 1975) 
d Radio flux at 1 GHz in the unit of 10~ 3 Jy. 

-l 



; Flux in the 2.0-10.0 keV band in the unit of 10 12 ergs cm 2 
Absorption column in the unit of 10 22 cm -2 . 
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Fig. 5. The XIS spectra of region 3 with the MEKAL model. 

Table 4. Best-fit parameters of the region 3 spectrum 



MEKAL 

Temperature [keV] 
Abundance a 
Normalization b [cm -3 ] 

phabs 
X 2 /d.o.f 

a Abundance ratio relative to the solar value (Anders & 
Grevesse, 1989). 

b Normalization is *° .> f n £ nijdV where Da is the 
angular size distance to the source (cm), n e and V are the 
electron density and the plasma volume. 
c Absorption column in the unit of 10 22 cm -2 . 

5. Timing Analysis 

In order to understand the nature of the point sources, we carried out timing analysis. 
5.1. Source A 

We searched the Chandra data for a pulsation from source A. Unfortunately, source A locates 
close to the chip boundaries of the ACIS-I, and is affected by the instrumental dithering effects. In 
fact, the source drops into the boundary every 1000 s, during which the source count diminishes 
significantly. Removing these time intervals with a Good Time Interval filtering, we made a light 
curve and search for a pulsation. The power spectrum on the basis of a 3.2 s binning light curve in 
the 0.3-10 keV band is shown in Fig. 6. We do not detect any pulsation in the period range between 
6.6-3000 s. 
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Fig. 6. Power spectrum of source A in the 0.3-10.0 keV band with the time bin size of 3.2 sec. 
5.2. Source B 

We created a light curve of source B in the 0.5-2.0 keV band with Suzaku, which is shown in 
Fig. 7. The time bin size is 512 s. Although we detected no drastic flare event, the light curve seems 



region3 source - back 512sec 



o 0.05 



f- 

< 



5*10 4 



10 5 
TIME s 



1.5*10 5 



Fig. 7. The light curve of region 3 in the 0.5-2.0 keV band. The time bin size is 512 s. 

to show flickering. In fact, Kolmogorov-Smirnov test indicates the probability of no variability in the 
light curve is 0.0012. 

6. Discussion 

6.1. Thermal Component 

We have calculated the electron number density and the age of the diffuse thermal plasma 
of region 1 and 2 on the basis of the best-fit parameters summarized in table 2 and 3. We assume 
that the plasma in region 1 distributes uniformly within a sphere with a radius of 1.'4, which is the 
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Half-Width-of-Half-Maximum (HWHM) radius obtained from the Chandra image (Aharonian et al. 
2008). Assuming the distance to CTB37B of 10.2 kpc (Caswell et al. 1975), we obtain the real radius 
of region 1 to be r rcg i — 1.3 x 10 49 cm. Accordingly, the volume of the region 1 plasma is V^ cg i — 
(4/3)7rr 3 cgl = 9.2 x 10 57 cm 3 . In the same way, from HWHM of an estimated image size of 1.'9 x 0'9, 
the real length of semi-major and semi-minor axes result in r reg2j i = 1.7 x 10 19 cm and r rcg 2 >s = 
8.2 x 10 18 cm, respectively. Assuming that the line-of-sight extent of region 2 is r reg 2, s , the volume of 
region 2 is calculated to be V rcg2 = (4/3)7rr re g2,irr eg 2 )S = 4.8 x 10 57 cm 3 (the resultant electron density 
becomes smaller by a factor of \/2 if we substitute r reg2) i for r rcg2iS ). With the aid of the emission 
measures obtained from the spectral fitting, E.M. = J n e nH<iV = 2.ll 1 ;o x 10 58 cm -3 (region 1) and 
2.6 x 10 57 cm -3 (region 2 upper limit, see § 4.2), the electron number density of region 1 and 2 are 

n eireg i = 1.7(1.2-2.2) [cm" 3 ] 

and 

™ e ,re g 2 < 0.82 [cm" 3 ], 

where we adopt the relation n e ~ 1.24 riu for fully ionized plasma. The parameter region in parenthe- 
sis is those allowed at the 90% confidence level. Assuming the strong shock, we obtain the pre-shock 
densities to be ~0.43 cm 3 and ^0.21 cm 3 , respectively, which are significantly lower than the aver- 
age density of the interstellar matter in the galactic plane. That means CTB37B exploded in a low 
density space. The density of region 1 together with the ionization parameter obtained from the fit of 
the region 1 spectra n e!Teg it = 3.5^{ 3 x 10 10 [cm" 3 s], enables us to estimate the age of the plasma 
observationally for the first time as 

t T egi = 6.5 (3.7-31) x 10 2 [yr]. 

CTB37B is one of the best candidates of SN393 in Chinese historical record (Stephenson & Green 
2002). The plasma age calculated from the observed ionization parameter and emission measure 
supports this identification. 

The number of electrons (N e = n e V) in region 1 and 2 are iV rcgl = 1.6 (1.1 — 2.0) x 10 58 and 
N e>Teg 2 < 3.9 x 10 57 . As a result, the total mass included in the two regions are 15 (11-2O)M and 
<3.7M , respectively, and the thermal energy (E = ^(N c + N H + N Uc )kT) are 

E vcgl = 6.4 (4.1 - 8.6) x 10 49 [ergs] 
E TCg2 < 1.6 x 10 49 [ergs]. 

under the assumption of energy equipartition between electrons and ions. The total thermal energy 
could be larger if other portions of the remnant are included and if the proton temperature is signif- 
icantly larger than the electron temperature as is expected for supernova remnants with large shock 
velocities (Ghavamian et al. 2007). 
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6.2. The Nature of the Point Sources 

6.2.1. Source A 

The best-fit spectral parameters of Chandra source A and those of Suzaku region 1 (source 
A and the diffuse thermal emission) are summarized in table 2. The measured hydrogen column 
densities (Nu - 4 x 10 22 cm -2 ) are the same between Chandra and Suzaku. However, since the 
Suzaku spectra below ~1 keV is dominated by the diffuse thermal component (Fig. 3), Nu obtained 
with Suzaku is determined mainly by the diffuse emission. This implies that source A is probably 
associated physically with the diffuse thermal emission in region 1, and is most likely a neutron star 
(or a black hole) born with the supernova explosion leading to CTB37B. In fact, the large photon index 
of 3.2 and luminosity of 2.2x 10 34 ergs s -1 are both consistent with those of anomalous X-ray pulsar 
(AXP; Fahlman & Gregory 1981; Kuiper et al. 2006). Young age of ~700 yr (§ 6.1) also supports 
this interpretation (e.g. the AXP 1E1841-045 is associated with Kes 73 whose age is estimated to be 
500-1000 yr; Tian & Leahy 2008). Since the spin period of AXP is in the range 6-12 s, it is natural 
that we have found no evidence of pulsation from the Chandra data, because the frame time of the 
ACIS-I is 2.3 s. In addition, the flux of the power law measured by Suzaku is apparently greater 
that that with Chandra by a factor of 1.8 (§4.1, table 2). Note that the point spread function of the 
Suzaku XRTs is not so sharp. Hence one may doubt that part of the power-law flux detected by 
Suzaku can be attributed to a putative diffuse non-thermal emission. However, ~80% of the Chandra 
region 1 photons above 3 keV, where the power-law component is dominant in the Suzaku spectra (see 
Fig. 3), originates from source A (table 1). Since the diffuse emission occupies only a small fraction 
in region 1, the flux difference between Chandra and Suzaku cannot be explained unless source A 
has really varied. This kind of long term variation is detected also from some other AXPs (Gavriil & 
Kaspi 2002; Kaspi et al. 2003). We therefore conclude based on the Chandra and Suzaku observations 
that source A is probably a new AXP. We need a fast timing observation with an imaging detector 
with a time resolution of less than 1 s to confirm our conclusion. 

6.2.2. Source B 

From the spectrum fitting, The hydrogen column density obtained from the fit to the Suzaku 
region 3 spectra iVn < 4 x 10 20 cm~ 2 (table 4) is much smaller than that obtained from region 1 
and 2. This result indicates that source B is a foreground source. The best-fit plasma temperature 
of kT ~ 0.5 keV is reminiscent of an active star. The existence of flickering (§5.2) supports this 
suggestion. 

6.3. Non-Thermal Component 

X-ray emission from CTB37B is composed of the diffuse thermal component (region 1) and 
the non-thermal component (region 2) as well as a point source (source A), as demonstrated in §3 
and §4. Hence, CTB37B now is the third SNR after RCW86 and Cas A that possesses the thermal 
and non-thermal X-ray emissions and TeV 7-ray emission all together. In addition, CTB37B is now 
the fifth SNR, following RCW86, Cas A, RX J1713.7-3946 and Vela Jr., from which non-thermal 
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Table 5. Comparison of non-thermal component with other SNRs. 



Target Name 


L x * LtbV ' 


LrevlLx 




rYev 


References + 


CTB37B (region 2) 


0.97 0.59 § 


0.61 


1.5 


2.3 


(1) 


RCW86 


3.8 0.55 


0.14 


3.1 


2.5 


(2)(3) 


Cas A 


110 0.21 


0.0019 


3.1 


2.4 


(4)(5) 


RXJ1713 


6.5 0.42 


0.06 


2.4 


2.2 


(6)(7) 


VELA Jr. 


-1.5xl0- 2 0.033 


-2.2 


2.6 


2.1 


(8) 


SN1006 


2.1 <0.15 


<0.1 


2.7 




(9)(10) 


* Unabsorbed flux in the 2- 

* Unabsorbed flux in the 1- 


-10 keV band in the unit of 10 34 
-10 TeV band in the unit of 10 34 


ergs s — L . 
ergs s — 1 . 









* (1) Aharonian et al. (2006); (2) Bamba et al. (2000); (3) Hoppe et al. (2007); (4) Helder et al. (2008); 
(5) Albert et al. (2007); (6) Slane et al. (1999); (7) Aharonian et al. (2004); (8) Aharonian et al. (2005); 
(9) Bamba et al. (2008); (10) Aharonian et al. (2005b) 

§ TeV 7-ray flux contributed to region 2 was calculated as a quarter of the total flux. 

radiation is detected both in X-ray and TeV 7-ray bands. The fluxes of the non-thermal emission of 
these five non-thermal SNRs are compared in table 5. The non-thermal diffuse component detected 
from region 2 has a remarkably flat X-ray spectrum with a photon index of 1.5. Since this photon 
index is equal to the typical radio photon index (energy index a = 0.5), the non-thermal emission 
spectrum can be considered as extending from the radio band smoothly to the X-ray band, thereby the 
roll-off energy results in as high as >15 keV (table 3). This roll-off energy is higher than any other 
SNR that is accompanied by the non-thermal X-ray and TeV 7-ray emission, such as <9 keV for 
RX J1713.7-3946 (Takahashi et al. 2008), ~0.23 keV for SN1006 (Bamba et al. 2008), ~0.87 keV 
for RCW 86 (Bamba et al. 2005). This indicates high electron acceleration efficiency in region 2. In 
addition, the density around region 2 is considered as being lower than in region 1, given that only the 
upper limit of the thermal emission is obtained (§ 6.1). The lower density may indicate higher shock 
velocity due to the smaller deceleration, which is consistent with the fact that the roll-off frequency is 
proportional to the square of the shock velocity (Aharonian & Atoyan 1999). 

It should be a matter of debate whether the TeV 7-ray emission and the non-thermal X-ray 
emission from region 2 are produced by the same population of electrons, since the images of these 
two bands shown in Fig. 2 are far from similar at first sight. The TeV 7-ray image is compatible with 
a shell with a radius of ~4'-6' due to limited spatial resolution of H.E.S.S. (Aharonian et al. 2008), 
which is compatible with the size of the radio shell. We thus assume that the TeV 7-ray emission is 
powered through 1-Zone Inverse Compton scattering (IC) of the cosmic microwave background due 
to the accelerated electrons. 

The maximum electron energy (E max ) can be evaluated by the shape of TeV 7-ray spectrum. 
Using the H.E.S.S. spectrum whose photon index is 2.3 (Aharonian et al. 2006) shown in Fig. 8 in 
blue, we determined -E max of 10 TeV. The red line in Fig. 8, on the other hand, is the X-ray power-law 
spectrum of region 2 with a photon index of 1.5. A series of the dashed plots are the model spectra 
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calculated under the assumptions of E max of 10 TeV, an index of the electron energy distribution of 
2, and various magnetic field (0.1, 1.0 and 10.0 fiG). It is clear from this figure that 1-zone IC model 
is unable to explain the synchrotron X-ray spectrum with any magnetic field strength. This suggests 
that TeV 7-ray emission is due to multi-zone IC scattering, or the decay of neutral pions generated by 
the high energy proton impacts. 
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Fig. 8. Spectrum energy distribution of region 2 from the X-ray and TeV 7-ray bands. The dashed plots on the 
left are the synchrotron radiation models with various magnetic field, and on the right is IC model spectrum. All 
the curves assume the maximum energy of electrons of 10 TeV. Blue line is the spectrum from H.E.S.S observation 
(T = 2.3; Aharonianet al. 2006), and the red line is the Suzaku region 2 spectrum which is apowerlaw with T = 1.5. 



7. Conclusion 

We have obtained with Suzaku the images and the high quality spectra of the supernova rem- 
nant CTB37B. The X-ray diffuse emission region coincides with that of radio and TeV 7-ray. The 
X-ray emission consists of thermal and non-thermal diffuse components as well as a point source re- 
solved by Chandra. CTB37B is the third SNR from which thermal and non-thermal X-ray emissions 
as well as TeV 7-ray emission are detected all together, and the fifth SNR that is accompanied by 
non-thermal emission both in X-ray and TeV 7-ray bands. 

The diffuse thermal emission can be best described by a non-equilibrium collisional ioniza- 
tion plasma model (NEI model) with a temperature, an ionization parameter (n e t [cm~ 3 s]), and the 
abundances of 0.9±0.2 keV, 3.5 t\\ x 10 10 , and ~ 0.5Z & (Mg, Si), respectively. The image size 
and the observed emission measure provides the number density of the thermal electrons before the 
shock to be 0.2-0.4 cm -3 , which is significantly lower than that of the Galactic plane. This suggests 
that the supernova explosion associated with CTB37B took place at a low density space. From the 
ionization parameter and the number density of the thermal electron, the age of the plasma is found to 
be ~650 till Y r - This is consistent with the tentative identification of CTB37B with SN393 within 
the error. 
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In contrast, the diffuse component occupying southern part of CTB37B (region 2) is non- 
thermal and represented by a power-law model or a srcut model. The photon index of 1 .5 is signif- 
icantly smaller than any other non-thermal SNR, but is consistent with that of a typical non-thermal 
SNR in the radio band. The srcut model fit with its normalization set free to vary therefore results in a 
high roll-off energy of >15 keV. Under the assumption that TeV 7-ray was emitted by 1-zone IC scat- 
tering, there are no solution for magnetic field strength that can reproduce the observed synchrotron 
spectrum in X-ray. This suggests that TeV 7-ray is produced by multi-zone IC scattering, or by the 
decay of neutral pions generated by the high energy proton impacts. 

Owing to the high spatial resolution of Chandra, a point source is resolved from the brightest 
part of the Suzaku image of CTB37B (region 1). Its association to the diffuse thermal emission 
indicated by A^h, the photon index of ~3, the X-ray luminosity of order 10 34 erg s _1 , and the long 
term flux variation evident from the Chandra and Suzaku observations all indicate that the point source 
is a new anomalous X-ray pulsar. A high speed photometric observation is encouraged. 

The authors are greatful to all of the Suzaku and H.E.S.S. team members. We also thank H. 
Yamaguchi, T. Tanaka and J. Vink for useful comments. This work was supported in part by Grant- 
in-Aid for Scientific Research of the Japanese Ministry of Education, Culture, Sports, Science and 
Technology, No. 19-4014 (A. B.), No. 18740153, No. 19047004 (R. Y). 
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